Evolution of second-order cosmological perturbations 
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We present a method for constructing gauge- invariant cosmological perturbations which are gauge- 
invariant up to second order. As an example we give the gauge-invariant definition of the second 
order curvature perturbation on uniform density hyp ersurf aces. Using only the energy conservation 
equation we show that this curvature perturbation is conserved at second order on large scales for 
adiabatic perturbations. 
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I. INTRODUCTION 

In this letter we report new results for the study 
of second-order perturbations about a Friedmann- 
Robertson- Walker (FRW) spacetime. The study of 
second-order perturbations to date has been limited by 
two main problems: firstly the difhculty in defining truly 
gauge-invariant perturbations at second order, and sec- 
ondly the complexity of the resulting Einstein equa- 
tions 0,0. 

We address the first of these challenges by extend- 
ing an approach to the construction of gauge-invariant 
variables previously advocated for first-order perturba- 
tions U]. We make an unambiguous physical definition 
of the perturbation, and by building this into the math- 
ematical description of the physical perturbation, con- 
struct a gauge-invariant quantity. As an example, we 
give a gauge-invariant definition of the curvature pertur- 
bation on uniform density hypersurfaces. 

We avoid much of the complexity of the second-order 
field equations by considering only the local energy con- 
servation equation 0, in the large-scale limit where 
we neglect all spatial derivatives. This gives us a simple 
result establishing the constancy of the large-scale curva- 
ture perturbation on uniform density hypersurfaces, up 
to and including second order, for adiabatic perturba- 
tions. 

To use, for example, the observed Gaussian dis- 
tribution of cosmic microwave background (CMB) 
anisotropics to test models for the origin of struc- 
ture in the very early universe, there is an implicit as- 
sumption that there is negligible growth of second order 
perturbations on large scales. Simple models of infiation 
produce an almost Gaussian distribution of density per- 
turbations J\ , and recent studies have shown that the co- 
moving curvature perturbation does not evolve at second 
order on large scales during slow-roll inflation '3, 0, ^3 ■ 
But non-linear evolution from the end of inflation up un- 
til the time of last-scattering of the CMB could produce 
deviations from a Gaussian distribution. The existence of 



a constant curvature perturbation at second order shows 
that the observed distribution of CMB anisotropics can 
be used to directly constrain the distribution of perturba- 
tions produced in the very early universe in simple infla- 
tion models that predict adiabatic density perturbations 
after inflation. 



II. SECOND-ORDER PERTURBATIONS 

Observations on scales close to the Hubble scale 
seem to be consistent with an almost homogeneous and 
isotropic spacetime that can be described by small per- 
turbations about a Friedmann-Robertson- Walker metric. 

Any tensorial quantities can be split into a homoge- 
neous background and inhomogcncous perturbation 

T(7y, x') = To(r/) + <5Ti(ry, x') + i<5T2(r,, x') + ... (2.1) 

where we use subscripts 1 and 2 to denote first and second 
order perturbations. 

We will consider perturbations about a spatially flat 
FRW background metric 

ds^ ^ c? {-d^^ + 8i^dx^dx^\ . (2.2) 

where r\ is conformal time and a = a(r\) is the scale factor. 
The metric tensor including second-order perturbations 
can be written as 

500 = -a' (1 + 201 <^2) , (2.3) 

ffo^ = a' [bu + \^B2^ , (2.4) 

= a" [(1 - 2^^! - V2) + 2Ciy + Cay] . (2.5) 

We will refer to '0 a-s the curvature perturbation as it 
describes the intrinsic scalar curvature of constant-?? hy- 
persurfaces on large scales. 

Perturbations can be split into scalar, vector, and ten- 
sor modes, according to their transformation behaviour 
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on spatial 3-hypersurfaces pjj . For instance we can write 



(2.6) 



where iJ is a scalar perturbation, Fi is a divergence-free 
vector, and hij a transverse, trace-free tensor perturba- 
tion. 



III. GAUGE TRANSFORMATIONS 

Under a second-order coordinate transformation 

^^ = ^^+ef + ^(C^i''+^2) ' (3-7) 

any tensor T and its perturbations defined in Eq. (|2.1|) 
transform as 0] 

6T2 = 5T2+L5,To + l2^To + 2%^Ti. (3.8) 

Thus under a first-order transformation — (ai, PI), 
a scalar quantity such as the density, p, transforms at 
first-order as 

Spi ^ 5pi + p'^ai , (3.9) 

while at second order, writing ^2 — (0^27/32); have 

Sp2 = Sp2 + p'oa2 + cti [pqUi + PqUi + 2Spi'] 

+ i2Sp,+p'^ai)Jl. (3.10) 

For the first order curvature perturbation we have 

^i^^Ji-Hai, (3.11) 

where Ti. = a' /a, while at second order we get from 
Eq. 

^2 = V2 - ai [Uai' + {n' + 2^2) ai - 2ij[ - AHt/ji] 
-Ha2 - (Hai - 21/^1 ),/31. (3.12) 



IV. GAUGE CHOICES AND 
GAUGE-INVARIANT VARIABLES 



In practice one can construct gauge-invariant defini- 
tions of unambiguous, that is physically defined, pertur- 
bations. These are not unique gauge-independent pertur- 
bations, but are gauge- invariant in the sense commonly 
used by cosmologists to define a physical perturbation. 
We draw a distinction here between quantities that are 
automatically gauge- independent, i.e., those that have 
no gauge dependence (such as perturbations about a 
constant scalar field), and quantities that are in gen- 
eral gauge-dependent (such as the curvature perturba- 
tion) but can have a gauge-invariant definition once their 
gauge-dependence is fixed (such as the curvature pertur- 
bation on uniform-density hypersurfaces). Although this 
approach has been widely used, at least implicitly, to con- 
struct gauge-invariant quantities at first-order P, 'l^ , it 
has not previously been used at higher-order. In this let- 
ter we show that it is possible to define gauge-invariant 
quantities at second-order corresponding to physical per- 
turbations. 



1. Uniform density hypersurfaces 

The uniform density hypersurfaces are defined by set- 
ting = to the required order in perturbation theory. 
We find that for a specific spatial gauge this leads to a 
specific temporal gauge to the required order. 

From Eq. (|3.9|) we see that setting Spi = to first- 
order requires a gauge shift from an arbitrary gauge 



Spi 
P'o 



(4.13) 



Leaving for the moment the spatial gauge dependence, 
we see from Eq. (|3.1U|I . and using Eqs. H4.13|l to fix the 
first-order temporal gauge shift, that to second order we 
require 



a2 



1 

To 



Sp2 - -^Spi'dpi + Spi 
Po 



(4.14) 



To completely fix the second-order temporal gauge 
shift H4.14|l picking out uniform density hypersurfaces we 
must also specify the first-order spatial gauge shift f]\. 
For example, a natural choice is to pick worldlines co- 
moving with the fluid. The fluid 3-velocity transforms as 



A gauge-invariant theory of linear perturbations about 
FRW metric was proposed by Bardeen and sub- 
sequentl y d eveloped by many authors (for example see 
Refs.|il|llllllllli3)- No such gauge- invariant formal- 
ism has been developed for non-linear cosmological per- 
turbations ^] . According to the Stewart- Walker lemma 
[Tsl any truly gauge-independent perturbation must be 
constant in the background spacetime. This apparently 
limits ones ability to make a gauge-invariant study of 
quantities that evolve in the background spacetime, e.g., 
density perturbations in an expanding cosmology. 



v'^v'-f3". (4.15) 

Thus from an arbitrary spatial gauge we can transform to 
the comoving gauge by the spatial gauge transformation 




(4.16) 



In this case there is a constant of integration correspond- 
ing to the arbitrary choice of spatial coordinates labelling 
the worldlines on an initial time-slice. 
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We are now able to construct gauge-invariant defini- 
tions for any metric or matter perturbations for comoving 
observers on uniform density hypersurfaces. As an exam- 
ple we give the curvature perturbation. To first-order, 
using Eqs. 13.11| and (|4.13() . we recover the well-known 
expression |4Lll9l| 



-Ci = ^1 



5pi 



(4.17) 



For the second order curvature perturbation on uniform 
density hypersurfaces, along comoving worldlines, we use 
Eq. |TT^ with Eqs. (jlT^ . and (jlTB . to give 



C2 



■02 



Po Po Po 



5pi 



n_i Po nil 0'l_/2 

o 1 /t— ; — ri — In. 



(4.18) 



We can also give gauge-invariant definitions for scalar 
quantities on uniform density hypersurfaces with comov- 
ing worldlines. As an example we write the pressure per- 
turbation on uniform-density hypersurfaces at first order 



5Pi 



P' 

Po 



(4.19) 



We can identify this as the usual gauge-invariant defini- 
tion of the non-adiabatic part of the pressure perturba- 
tion. At second order we have 



where the extra term on the right-hand-side of the ex- 
pression for the second-order non-adiabatic pressure per- 
turbation arises from the local variation of the adiabatic 
sound speed, = dP/dp. Note that the non-adiabatic 
pressure perturbation is automatically gauge-invariant 
at first order, but the second order non-adiabatic per- 
turbation is only gauge-invariant if the first order non- 
adiabatic pressure perturbation vanishes 2]. 

Thus we can write the gauge-invariant pressure per- 
turbation on uniform density hypersurfaces as 



SPi\ 



SPi 



nad 



SP2\p = 5P2„ad - 2^(5Pi; 

Pa 



ad ' 



(4.22) 
(4.23) 



We see that the pressure perturbation will vanish on 
hypersurfaces of uniform density for adiabatic perturba- 
tions. 



2. Uniform curvature hypersurfaces 

Instead of defining quantities on uniform density hy- 
persurfaces we can choose to work with uniform curva- 
ture slices. In some scenarios this can have the advantage 
of staying non-singular even when the uniform density 
hypersurfaces become ill-defined [2(iLl2lj |. 

We can define uniform curvature hypersurfaces by 
"01 = and -02 = 0, which fixes the temporal gauge 
shift and Ei = Q and Fl = to fix the spatial gauge shift 
to first-order. This implies for the first order temporal 
gauge shift, using Eq. (|3.11|) . 



"1 = TTT 



01 
H 



(4.24) 



6P2 



SP2 ~ ^Sp2 + 2 (sPi - ^Spi 
Po \ Po 



-Po 



5 Pi SPi'\ 26 p 



P'o 



P' 

^0 



P'o 



P" 
zJL 

P' 
^0 



Pi 
p'o 



5pI_ 
P'S 
(4 



and at second order, using Eq. H3.12|l . 



a2 



n 



n 



(4.25) 



These results are readily extended to systems involving 
scalar fields, as scalar fields obey the same transforma- 
tion rules as the energy density, given in Eqs. H3.9|) and 
(|3.1U|I . Hence one can write down the comoving cur- 
vature perturbation (i.e., the curvature perturbation on 
uniform scalar field hypersurfaces) or the relative entropy 
perturbation between two fields. 

For adiabatic perturbations the local pressure is a 
unique function of the local density, P = P{p)- Hence 
we can identify the non-adiabatic part of the pressure 
perturbation, to first and second orders, as 



20)vhere we used Eq. (|4.24|) and fix the spatial gauge 
shift 



<5-Plnad = SPi - clSpi , 
<5-P2nad = SP2 - clSp2 - 



dpo 



6pi 



(4.21) 



0{ = E,; 



(4.26) 



The density perturbation on uniform curvature hyper- 
surfaces is then, at first order, using Eqs. (|3.9|) . 



Spi 



Spi + ^A. 



(4.27) 



At first order the curvature perturbation defined on uni- 
form density hypersurfaces and the density perturbation 
on uniform curvature hypersurfaces are simply related by 



6pi 



Po 
H 



^1 



(4.28) 
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From Eq. (|3.10|) we get the definition of the second 
order density perturbation on uniform curvature hyper- 
surfaces 



n n 



p'on'\ 2 



-2^v^;^i + ^^iSpi' 



(4.29) 



Again, these results are readily extended to scalar fields 
which obey the same gauge transformations rules as the 
energy density. 



V. EVOLUTION OF CURVATURE ON LARGE 
SCALES 

Having defined gauge-invariant second-order variables, 
we now turn to finding evolution equations for these 
quantities. 

A fundamental question in cosmology is how density 
perturbations evolve in the large-scale regime where grav- 
itational perturbations cannot be neglected. In particu- 
lar it is important to establish whether non-linear evolu- 
tion could introduce significant non-Gaussianity, e.g., in 
CMB anisotropics, even in infiationary models where the 
large scale structure of the universe is supposed to arise 
from purely Gaussian fiuctuations at very early times. 

Despite the complexity of the field equations at second 
order (see e.g. Ref. (225) it is sufficient to use the local 
conservation of energy-momentum to establish the con- 
servation of if we neglect spatial gradients, which we 
expect to be valid on sufficiently large scales. This strat- 
egy was first employed in Ref. Q to establish the con- 
servation of C at first order. We will also neglect terms 
quadratic in first-order vector and tensor perturbations 
which can be calculated using the first-order equations of 
motion and can be shown to be small on large scales in an 
expanding universe. We defer a full treatment including 
the effect of these quadratic terms to future work. 

We define the energy momentum tensor as 

T^" = {p + P) u^u" + Pg"" + n"" , (5.30) 

where n'^" is the trace-free anisotropic stress tensor, and 
m'^ is the fiuid 4-velocity. 

Energy conservation is given by MjyV^T^" = 0. In the 
homogeneous background we have 



Po + 3-(po + ^o) = 0, 
a 

while dropping spatial gradient terms we obtain 



(5.31) 



Spi' + 3ni6pi+SPi)~3ipo + Po)ij[ ^ 0,(5.32) 
Sp2' + m {5p2 + 5P2) - 3 (po + Po) V^2 

~Q4,[[5pi+5Pi+2{pa + Po)i:i] ~ 0,(5.33) 



to first and second order respectively. These can be writ- 
ten in terms of the gauge-invariant curvature perturba- 
tion, ( defined in Eqs. (|4.17|) and (|4.18|l . giving 



and 



n 



H 



ip + P) 



SPi\ 



(5.34) 



(P + P) 
2 



Po 



SP2\p 

SPi\p-2{po + Po)Ci] Ci, (5.35) 



where the gauge-invariant pressure perturbation on uni- 
form density hypersurfaces is given in terms of the non- 
adiabatic pressure perturbation in Eqs. I|4.22|l and 1)4. 23|) . 
Thus the curvature perturbation on uniform density hy- 
persurfaces is constant at first and second order for adia- 
batic perturbations on large scales if we can neglect spa- 
tial gradients [25l |. 

Hence we expect an initially Gaussian distribution of 
adiabatic curvature perturbations in the very early uni- 
verse will remain Gaussian in the large scale limit. By 
the same token, primordial non-Gaussianity of the per- 
turbations may be indicative of non-adiabatic evolution 
in the early universe. 



VI. CONCLUSION 

In summary, we have given a procedure for defining 
gauge-invariant cosmological perturbations at first and 
second order. As an example we have given a gauge- 
invariant definition of the curvature perturbation on uni- 
form density hypersurfaces. We expect that this prescrip- 
tion could be extended to higher orders if desired. 

We were then able to show, using only the local energy 
conservation equation, that the curvature perturbation 
remains constant on large scales for adiabatic perturba- 
tions where we neglect spatial gradients. As shown in 
Ref. , a conserved perturbation exists of any quantity 
that obeys an autonomous local conservation equation. 
Thus we can construct conserved perturbations of the lo- 
cal energy density of any fluid with a barotropic equation 
of state. 

Note added: While writing up this letter, a non- linear 
result for the constancy of C in a long- wavelength approx- 
imation jj was reported by Rigopoulos and Shellard [2^ . 



Acknowledgments 

The authors are grateful to Marco Bruni, Kouji Naka- 
mura, David Lyth, David Matravers and Toni Riotto 
for useful discussions. This work was supported by 
PPARC grant PPA/G/S/2000/00115. KM is supported 
by a Marie Curie Fellowship under the contract number 
HPMF-CT-2000-00981. DW is supported by the Royal 



5 



Society. Algebraic computations of tensor components were performed using the GRTensorll package for Maple. 



[1] D. S. Salopek and J. R. Bond, Phys. Rev. D 42 (1990) 
3936. 

[2] M. Bruni, S. Matarrese, S. MoUerach and 
S. Sonego, Class. Quant. Grav. 14 (1997) 2585 



|rarXiv:gr-qc/9609040 
and M. Bruni, Phys. 



S. Matarrese, S. MoUerach 
Rev. D 58, 043504 (1998) 
arXiv:astro-ph/9707278 ; M. Bruni, F. C. Mena and 
R. K. Tavakol, Class. Quant. Grav. 19, L23 (2002) 
|arXiv :gr-qc/0107069 . 
[3] iTX Mahk and D. Wands, |arXiv:gr-qc/9804046| 
[4] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, 
Phys. Rev. D 62 (2000) 043527 arXiv:astro-ph/0003278 . 
[5] D. H. Lyth and D. Wands, Phys. Rev. D 68 (2003) 

103515, arXiv:astro-ph/0306498 . 
[6] E. Komatsu et ai, arXiv:astro-ph/0302223 
[7] A. Gangui, F. Lucchin, S. Matarrese and 

MoUerach, Astrop hys. J. 430, 447 (1994) 

arXiv:astro-ph/9312033 ; A. Gangui, Phys. Rev. D 
50, 3684 (1994) arXiv:astro-ph/9406014 ; L. M. Wang 
and M. Kamionkowski, Phys. Rev. D 61, 063504 (2000) 
arXiv:astro-ph/9907431 . 
[8] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, 
Nucl. Phys. B 667 (2003) 119, arXiv:astro-ph/0209156 . 
[9] J. Maldacena, JHEP 0305 (2003) 013 

[10] 

[11] 
[12] 

[13] 
[14] 

[15] 



arXiv:astro-ph/0210603 




G. Rigopoulos, 


arXiv:astro-ph/021214ll 



H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78 
(1984) 1. 

K. A. Malik, arXiv:astro-ph/0101563 . 

V. F. Mukhanov, H. A. Feldman and R. H. Branden- 

berger, Phys. Rept. 215 (1992) 203. 

R. Durrer, Fundamentals of Cosmic Physics, 1994, Vol. 



15., No. 3, pp. 209-339, astro-ph/9311041 . 
[16] J. c. Hwang and H. Noh, Gen. Rel. Grav. 31, 1131 (1999) 

arXiv:astro-ph/9907063 . 
[17] K. Nakamura, Prog. Theor. Phys. 110 (2003) 723, 

arXiv:gr-qc/0303090 . 
[18] J. M. Stewart and M. Walker, Proc. R. Soc. A, 341 49, 

(1974). 

[19] J. M. Bardeen, DOE/ER/40423-01-C8 Lectures given at 
2nd Guo Shou-jing Summer School on Particle Physics 
and Cosmology, Nanjing, China, Jul 1988. 

[20] K. A. Malik, D. Wands and C. UngareUi, Phys. Rev. D 
67 (2003) 063516 arXiv:astro-ph/0211602 . 

[21] F. FineUi and R. H. Brandenberger, Phys. Rev. Lett. 82, 
1362 (1999) arXiv:hep-ph/9809490 . 

[22] H. Noh and J. Hwang, arXiv:astro-ph/0305123 

[23] G. I. Rigopoulos and E. P. Shellard, Phys. Rev. D 68 
(2003) 123518, arXiv:astro-ph/0306620 . 

[24] After completing this work, our attention was drawn to 
a recent paper by Nakamura [l^ where a general proce- 
dure is proposed for constructing gauge-invariant pertur- 
bations upto third order of a generic spacetime. 

[25] An alternative argument for the existence of a conserved 
quantity at second order was given recently in Ref. 
defined in terms of the density perturbation on uniform 
expansion hypersurfaces. One can verify that for adi- 
abatic perturbations the second-order part of the con- 
served quantity defined in Eq. (11) of Ref. coincides 
with C,2 — 2("i , using the gauge-invariant defintions of ("i 
and ^2 given in Eqs. 14.1711 and 14.1811 . written in terms 
of the density perturbation evaluated on spatially flat 
hypersurfaces. 



